home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / QROMB.DEM < prev    next >
Text File  |  1991-04-29  |  791b  |  41 lines

  1. PROGRAM d4r4(input,output);
  2. (* driver for routine QROMB *)
  3. CONST
  4.    pio2=1.5707963;
  5.    n=5;
  6. TYPE
  7.    glnarray = ARRAY [1..n] OF real;
  8. VAR
  9.    glit : integer;
  10.    a,b,s : real;
  11.  
  12. FUNCTION func(x: real): real;
  13. (* Test function *)
  14. BEGIN
  15.    func := sqr(x)*(sqr(x)-2.0)*sin(x)
  16. END;
  17.  
  18. FUNCTION fint(x: real): real;
  19. (* Integral of test function func *)
  20. BEGIN
  21.    fint := 4.0*x*(sqr(x)-7.0)*sin(x)
  22.       -(sqr(sqr(x))-14.0*sqr(x)+28.0)*cos(x)
  23. END;
  24.  
  25. (*$I MODFILE.PAS *)
  26. (*$I TRAPZD.PAS *)
  27.  
  28. (*$I POLINT.PAS *)
  29.  
  30. (*$I QROMB.PAS *)
  31.  
  32. BEGIN
  33.    a := 0.0;
  34.    b := pio2;
  35.    writeln ('Integral of func computed with QROMB');
  36.    writeln;
  37.    writeln ('Actual value of integral is',fint(b)-fint(a):12:6);
  38.    qromb(a,b,s);
  39.    writeln ('Result from routine QROMB is',s:12:6);
  40. END.
  41.